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Failed  Elements  in  a  Phased  Array  Antenna 
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Abstract — Maintaining  low  spatial  sidelobes,  despite  the  failure 
of  individual  antenna  elements,  is  important  for  many  radar  and 
communication  systems  utilizing  phased  arrays.  Techniques  in 
the  literature  require  accurate  antenna  element  pattern  data.  In 
this  paper  we  present  a  new  algorithm  for  computing  low- 
sidelohe  heamforming  which  only  requires  the  original 
heamforming  weights  which  produce  low  sidelohes  when  all 
elements  are  functioning  normally.  The  algorithm  does  not 
require  accurate  knowledge  of  the  antenna  element  patterns,  and 
permits  user  adjustment  of  the  trade-off  between  sidelobe  level, 
taper  loss,  and  mainbeam  width.  Near  optimum  low  sidelobes  are 
demonstrated  in  several  examples. 

Index  Terms  —  Array  signal  processing,  beams,  linear  algebra, 
phased  arrays,  shaped  beam  antennas. 

L  INTRODUCTION 

Eor  many  phased  array  antenna  applications,  low  spatial 
sidelobes  are  required,  and  it  is  desirable  to  maintain  these  low 
sidelobes  despite  the  failure  of  one  or  more  antenna  elements. 
Eor  most  systems  the  array  manifold  (i.e.,  the  individual 
complex  antenna  element  patterns)  are  known  well  enough  to 
compute  a  low  spatial  sidelobe  pattern  despite  the  failure  of 
several  elements,  as  described  in  [1-8].  Eor  this  paper  it  is 
assumed  the  low  sidelobe  beamformers  with  no  failed 
elements  are  known.  However  with  failed  elements  the  array 
manifold  is  not  known  accurately  enough  to  apply  any  of  the 
referenced  approaches.  The  complex  antenna  element  patterns 
are  different  and  the  measurement  errors  are  too  large. 
Nothing  has  been  found  in  the  literature  for  this  problem.  The 
authors  published  a  paper  in  ICASSP  2012  on  this  topic  [9]. 
This  expanded  paper  provides  detailed  mathematical 
justification,  performance  bounds,  and  additional  techniques 
not  included  in  the  earlier  paper. 

Specifically,  consider  an  n  element  phased  array  and  the 
following  length  n  vectors.  Let  0  be  an  angular  direction 
which  for  a  two  dimenstional  array  is  a  vector  of  length  2  and 
for  a  one  dimentional  array  is  a  scalar.  Let  0  be  the  desired 
direction  for  a  low  sidelobe  beamformer.  We  assume  with  no 

Distribution  A:  Public  Release.  This  work  is  sponsored  by  the  Department 
of  the  Air  Force  under  Air  Force  contract  #FA8721-05-C-0002.  Opinions, 
interpretations,  conclusions  and  recommendations  are  those  of  the  authors  and 
are  not  necessarily  endorsed  by  the  United  States  Government. 

*S.  1.  Krich  is  with  Massachusetts  Institute  of  Technology  Lincoln 
Laboratory,  Lexington,  MA  02421  USA  (e-mail:  sik23@cornell.edu). 

1.  Weiner  is  with  Harvard  University,  Cambridge,  MA  02138  USA  on 
leave  from  Lincoln  Laboratory,  Lexington,  MA  02421  USA  (e-mail: 
iweiner@seas.harvard.edu). 


failed  elements  the  following  length  n  vectors  are 
UNKNOWN  and  KNOWN  and  all  vectors  except  for  £(0)  are 
unit  normed. 

UNKNOWN  for  all  0 

•  Vt(0)  :  true  steering  vector  to  angle  0  (i.e.,  on  the  array 
manifold) 

•  f(0):  antenna  array  manifold  (calibration)  errors 

KNOWN 

•  w(0):  low  sidelobe  beamformers  for  all  0  within  a 
neighborhood  around  0,  where  the  inner  product 
{w(0),Vf.(0 A))  is  small  for  0  +  A  in  the  sidelobe 
region 

•  V„(0)  =  [vM  +  £(d)]/  \vM  +  £(0)|  :  assumed 
steering  vector  to  angle  0 

The  foundation  of  this  paper  is  the  observation  that  with 
failed  elements  and  perfect  knowledge  of  the  array  manifold, 
a  nearly  optimal  low  sidelobe  beamformer  is  represented  by  a 
linear  combination  of  low  sidelobe  beamformers  with  no 
failed  elements,  w(0)’s,  in  a  neighborhood  around  0  under 
the  constraint  that  the  linear  combination  has  zeros  at  the 
location  of  the  failed  elements.  Eurthermore  since  each  of 
these  low  sidelobe  beamformers  achieve  low  sidelobes,  we 
would  expect  that  linear  combinations  of  them  in  a 
neighborhood  around  0  would  also  have  low  sidelobes. 

The  algorithms  in  this  paper  exploit  this  observation  using 
only  the  KNOWN  vectors.  The  paper  derives  equations  for 
estimating  taper  loss,  sidelobe  level,  and  mainbeam  region 
width,  and  the  algorithm  permits  optimization  trade-offs 
among  these  three  important  criteria.  Examples  achieve  a  near 
optimal  solution  comparable  to  accurately  knowing  v^(9)  for 
all  0.  Thus  for  some  systems,  in  particular  where  a  rapid 
recalculation  of  the  beamformer  with  failed  elements  is 
required,  it  may  be  preferable  to  use  this  method  even  if  Vt(0) 
is  accurately  known. 

The  reader  may  be  wondering  how  the  w(0)’s  could  have 
been  determined  in  the  first  place  without  knowing  the  array 
manifold.  Examples  are  airborne  radars  that  during  low 
sidelobe  antenna  pattern  calibration  flights,  form  beams  by 
minimizing  the  sidelobe  ground  clutter.  The  method  in  [10] 
iteratively  adjusts  the  beamformer  and  only  requires 
observability  of  the  single  beam.  Space  Time  Adaptive 
Processing  (STAP)  [11]  requires  a  fully  digitized  array.  STAP 
could  also  be  used  with  failed  elements  but  may  not  be 
available  during  normal  radar  operation. 
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Depending  upon  the  nature  of  the  antenna  element  failure, 
the  mutual  coupling  between  elements  can  change,  which 
impacts  the  true  array  manifold.  The  results  here  and  [1-8] 
apply  to  cases  in  which  these  changes  are  small.  Eor  example, 
failures  that  do  not  change  the  impedance  at  the  antenna  port 
have  a  negligible  effect  upon  the  mutual  coupling. 

Section  II  describes  methods  for  recomputing  the 
beamformer  weights  with  failed  elements.  Section  III  provides 
examples  for  both  one-  and  two-dimensional  arrays.  Section 

IV  applies  these  techniques  to  arrays  with  no  failed  elements 
to  modify  either  the  taper  loss  or  the  sidelobe  level.  It  is 
connnon  for  techniques  developed  for  phased  antenna  arrays 
to  also  be  applicable  in  the  temporal  domain  and  thus  Section 

V  discusses  how  the  same  technique  is  applicable  for  pulse- 
Doppler  radars  with  missing  pulses  or  missing  range  samples 
due  to  interference.  Section  VI  is  a  sunnnary. 

II.  Tailed  element  algorithm 

Eig.  1  shows  a  digitally  controlled  beamformer,  w(0), 
applied  to  the  n  element  array  to  produce  a  single  beam.  If 
elements  are  digitized  the  beamformer  is  digital;  otherwise  it 
is  analog. 

A.  Beamformer  with  good  antenna  calibration 

We  start  by  describing  a  procedure  which  could  be  used 
with  good  knowledge  of  the  array  manifold.  With  no  failed 
elements,  we  could  compute  a  low  sidelobe  beamformer  w(0) 
by  modeling  a  covariance  matrix  with  sidelobe  interference. 

M>(0)  =  n[Yl  +  (1  -  Y)M(e)WVa(d) 

~  1  (^) 

M(0)  =l  I,  Va(et)Va(etr  . 

OisasL 

where  I  is  an  identiy  matrix  representing  the  thermal  noise, 
M(0)  is  the  modeled  sidelobe  interference  covariance  matrix 
with  no  noise,  y  controls  the  mixture  of  modeled  interference 
to  thermal  noise,  is  the  set  of  sidelobe  directions,  L  is 
the  number  of  terms  in  the  sum,  is  a  normalizing  scale 
factor,  and  H  is  Hermitian  transpose.  With  failed  elements, 
the  corresponding  rows  and  columns  of  vectors  and  matrices 
are  deleted.  With  very  good  knowledge  of  the  array  manifold, 
this  procedure  works  very  well  but  unfortunately  with  or 
without  failed  elements,  the  procedure  will  fail  to  achieve  low 
sidelobes  if  the  array  manifold  errors  are  too  large. 

Using  just  the  KNOWN  quantities  in  Section  I,  the  next 
three  sections  derive  procedures  for  selecting  the  number  of 
w(0)’s  to  include  in  this  linear  combination,  their  relative 
weighting,  plus  estimates  for  the  resulting  sidelobe  level 
achieved  and  the  taper  loss.  The  ability  to  estimate  sidelobe 
level  and  taper  loss  permits  optimization  trades. 

B.  Beamformer  with  constraints 

Here  we  modify  (1)  by  constraining  the  solution  to  a  specific 
subspace.  Let  M4:  be  a  matrix  whose  columns  are  K  low 
sidelobe  beamformers  surrounding  the  direction  0  of  interest. 
These  should  be  closely  spaced  (less  than  a  beamwidth)  and 
preferably  include  w(0)  as  a  column.  We  constrain  the 


solution  to  the  vector  space  spanned  by  M4:  as  follows.  Let 
J  <  K  be  the  number  of  failed  elements  and  D  = 

[d^,d2, ...  ,djY  be  a  vector  for  the  locations  of  the  failed 
elements.  Within  the  space  spanned  by  M4:  is  a  subspace  S  of 
dimension  K  —  ]  where  all  vectors  in  S  have  a  0  at  the 
locations  of  the  failed  elements.  We  compute  this  subspace  as 
follows.  Using  MATLAB  notation,  is  3.  J  X  K 

matrix  of  only  the  rows  of  with  failed  elements.  The 
K  X  K  —  J  matrix 

Z  =  null(VK^(D,;))  (2) 

is  an  orthonormal  basis  for  the  null  space  of 
obtained  from  the  singular  value  decomposition.  That  is, 
Wf^(D,:)Z  is  a  J  X  (K  —  J)  matrix  of  Os  and  thus  M4:Z 
is  an  n  X  (K  —  J)  matrix  with  Os  along  the  rows 
corresponding  to  the  location  of  the  J  failed  elements.  We 
constrain  the  solution  to  subspace  S  thereby  modifying  (1)  as 
follows 

w{e)=fiw^z[(w^zr[ri  +  (i  -  r)M{e)](w^z)Y\w^zrva(0) 
=fiW^z[Yzxw«w^z  +  (1  -r:izXG(e)zY\w„zrva(0) 

where 

G(e)  =  w^M(e)WK.  (4) 

The  next  two  sections  develop  different  approximations  to 
G(^0)  based  upon  different  assumptions  on  how  the  low 
sidelobe  beamformers,  w(0)’s,  were  formed. 

C.  Uncorrelated  beam  approximation  to  G 

We  never  actually  compute  M(0)  but  it  is  helpful  to  think  of 

D.SL  the  sidelobe  region  which  is  outside  the  mainlobe  of  all 
the  beamformers  in  Wf^,  Erom  (1)  and  (4)  we  can  write 

G(o)  =  j  y  [wi^Va(di)][Va(etrw^]  .  (5) 

OisasL 

where  is  a  vector  for  the  sidelobe  level  at  9i  for 

each  of  the  K  beamformers  in  W^.  We  assume  now  that  the 
process  by  which  all  of  the  beamformers  in  M4:  were  produced 
yield  the  same  average  sidelobe  level  a,  but  the  actual 
sidelobes  averaged  over  flsi  approximately  uncorrelated 
from  beam  to  beam.  Eor  example  if  w(0)  =  w(0)  +  x(0) 
where  w(0)  has  lower  sidelobes  than  w(0)  and  the  sidelobe 
level  is  established  by  random  vector  x(0)  which  is  different 
for  each  0.  Thus  due  to  the  averaging  over  the  sidelobe  angles, 
the  off  diagonal  terms  in  (5)  are  small  and  can  be 
approximated  as  zero.  Our  assumption  yields  the 
approximation 

G^(G)  =  al  .  (6) 

D.  Covariance  estimate  to  G 

Any  small  angle  dependent  variations  in  the  element  antenna 
patterns  from  element  to  element  will  impose  a  fundamental 
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minimum  on  the  level  of  antenna  sidelobes  that  can  be 
achieved.  Here  we  consider  the  case  where  the  w(0)’s  are 
close  to  this  limit  with  small  taper  loss  as,  for  example,  if  they 
were  generated  using  (1)  with  y  =  0  with  excellent  knowledge 
of  the  array  manifold. 

Under  the  conditions  described,  the  original  beamformer 
w(0)  is  close  to  optimal  with  respect  to  G(0),  so  it  is 
reasonable  to  seek  an  approximation  62(0)  which  preserves 
this  property;  This  can  be  expressed  via  covariance  inversion, 
or  by  the  equivalent  requirement  that  ^2(0)^  (0)  OC  Vaid).  In 
addition,  similarly  to  above,  we  require  that  sidelobe 
estimates  for  each  of  the  Wj^  agree  with  observation  [see  (9) 
below]. 

A  straightforward  computation  shows  that  the  following 
(positive-definite,  Hermitian)  matrix  satisfies  these  two  key 
model  requirements: 


G2(0)  =  a 


Va(0)  W 


■w^vaie)vaiefw„ 


+  diag  1 


(7) 


where  Wj^,  k  =  1,2,  ...K  are  the  chosen  beamformers  which 
form  the  columns  of  1/14:,  the  MATLAB  notation  for  a 
diagonal  matrix  is  being  used. 

E.  Taper  loss  and  sidelobe  estimate 

In  choosing  parameters  K  and  7  it  is  important  to  have  an 
estimate  for  the  taper  loss  and  changes  in  the  average 
sidelobes.  The  taper  loss  estimate  is 

TL(wie))  =  <  1  (8) 

where  w(0)  and  Va(0)  are  unit  normed.  The  manifold  errors 
may  be  too  large  to  achieve  the  desired  sidelobes  using  the 
methods  in  [1-8];  however,  this  taper  loss  estimate  is  quite 
accurate  unless  the  manifold  errors  are  very  large. 

We  can  write  w(0)  =  Wj^c  where  c  is  determined  from  (3). 
Since  M(0)  is  a  covariance  matrix  of  steering  vectors  in  the 
sidelobe  region,  the  average  sidelobes  is 

SL{w(d))  =  w(e)^  M(e)w(d) 

=  c^W^M(e)WKC  (9) 

=  c^G(^0)c  . 

Using  the  approximation  Gi(^0)  in  (6)  or  (7)  in  place  of  G(^0), 
we  estimate  the  change  of  sidelobes  as 

~  .  c^Gi(d)c  .-1^. 

A5Lest(iv(0))  = - - -  i  =  1  or  2  .  (10) 

One  does  not  need  to  know  the  value  for  a  to  use  (10)  since  by 
the  way  Gi(^0)  is  defined  in  (6)  or  (7)  a  cancels.  The  appendix 
develops  some  rough  upper  bounds  on  the  error  incurred  by 


using  Gi(^0)  .  Eor  small  enough  perturbations  from  the  original 
vector,  (9)  gives  a  good  approximation  even  if  (6)  or  (7)  are 
not  uniformly  good  approximations  to  G(^0). 

F.  Choosing  parameters  K  and  y 

In  (3)  7  controls  the  mixture  of  modeled  interference  to 
thermal  noise.  Eor  7  =  1  this  is  a  projection  of  Va(0)  onto 
the  space  spanned  by  the  columns  of  Wj^Z  and  thus  has  the 
best  taper  loss  but  the  highest  estimated  sidelobes.  As  y 
decreases,  the  taper  loss  monotonically  degrades  and  the 
estimated  sidelobe  power  monotonically  improves.  (See,  for 
example,  theorem  II  in  [12].)  7  <  1  has  the  effect  of 

regularizing  the  matrix  Z^W/^Wj^Z  by  reducing  the 
contribution  of  the  eigenvectors  corresponding  to  its  small 
eigenvalues.  These  monotonicity  properties  are  very  useful  for 
choosing  parameters  in  the  optimization  procedures  below. 

As  is  always  the  case  when  choosing  beamforming  weights, 
trade-offs  can  be  made  between  taper  loss,  sidelobe  level,  and 
mainbeam  region.  Therefore,  parameter  choices  must  be 
determined  in  part  by  the  properties  most  important  to  the 
application.  However,  the  choice  is  simplified  greatly  by  the 
monotonicity  observations  above.  K  directly  impacts  the  width 
of  the  mainbeam  region.  Under  our  correction  scheme,  the 
coefficients  of  the  linear  combination  of  neighboring  beams 
are  roughly  the  corrected  pattern  gain  at  their  corresponding 
look  directions.  Therefore  a  larger  K  will  create  a  wider 
mainbeam  region  manifested  as  some  combination  of  a  wider 
mainlobe  and/or  high  first  sidelobes.  Note  that,  even  with 
antenna  manifold  errors,  the  shape  of  this  resultant  mainbeam 
region  can  be  predicted  and  possibly  adjusted  depending  on 
the  needs  of  the  application.  No  matter  which  method  is 
chosen  to  estimate  G(^0),  two  procedures  follow  depending 
upon  the  desired  goals.  Eor  both  procedures  it  is  convenient  to 
take  a  =  1  since  we  only  estimate  the  change  in  the  sidelobe 
levels  in  (10)  and  the  actual  sidelobe  levels  may  be  unknown. 

In  Procedure  1  we  set  a  goal  for  the  change  in  the  average 
sidelobes,  to  achieve  and  a  bound  on  the  taper  loss, 

TLbound  to  avoid  exceeding.  These  are  estimated  using  (8)  and 
(10).  It  is  usually  impossible  to  achieve  both  of  these  exactly. 
Procedure  2  is  the  reverse  with  a  bound  on  the  change  in 
average  sidelobes,  and  a  goal  for  the  taper  loss, 

^^goal- 

Procedure  1  is  flow  charted  in  Pig.  2  in  which  values  for  K 
and  7  are  found  such  that  ASL^sti}^)  —  ^^^goai  the  taper 
loss,  TL(w),  is  lower  bounded  by  TLbound-  (Note:  ASlg^^^i  =  1 
means  unchanged  sidelobes  and  TL(w)  is  negative.)  Pinding  7 
in  the  interval  [0,1]  is  simplified  by  the  monotonicity 
discussed  above  for  the  function  in  (8).  K  is  increased  until  an 
acceptable  criteria  are  met. 

Procedure  2  is  flow  charted  in  Pig.  3  in  which  taper  loss 
TL(w)  =  TLgoai  and  ASL^stC^)  ^  bound-  Here  again  K  is 
increased  until  an  acceptable  criteria  are  met. 

With  either  procedure,  should  the  goal  and  bound 
requirements  be  too  stringent  to  achieve,  K  will  continue  to 
increase  to  an  unacceptably  large  value  and  a  solution  may  not 
be  found,  in  which  case  the  goal  and  bound  should  be  reset 
and  the  procedure  repeated.  K  is  constrained  to  be  odd  to 
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provide  symmetric  beams  around  the  beam  of  interest  though 
an  odd  value  is  not  an  inherent  requirement. 

III.  Examples 

A.  64  element  uniform  linear  array  with  half  wavelength 

spacing 

The  authors  feel  it  is  important  that  the  validity  of 
algorithms  such  as  the  one  described  in  this  paper  be  tested 
under  real-world  conditions.  Specifically,  we  test  the 
procedure  using  a  perfect  uniform  linear  array  as  the  assumed 
array  manifold,  but  a  true  array  manifold,  Vt(0), 

having  small  perturbations  from  the  perfect  uniform  linear 
array.  The  manifold  errors,  f(0),  limit  the  achievable 
beamformer  sidelobes  using  (1)  to  -30  dB.  However,  the 
beamforming  weight  vectors,  w(0),  to  achieve  low  sidelobes 
are  available.  The  beams  in  Wj^  are  spaced  by  a  half 
beamwidth  (i.e.,  sin(0j)  =  i/n). 

The  manifold  errors,  f(0),  are  modeled  as  independent 
circularly  Gaussian  on  each  element  of  the  array.  Since  on 
individual  elements  the  errors  are  likely  to  be  correlated  with 
look  direction,  for  a  beamwidth  change  in  look  direction  the 
errors  are  correlated  by  0.98  on  each  element  but  independent 
from  element  to  element.  This  is  achieved  by  starting  with  an 
nxm  matrix  of  independent  circular  Gaussian  random 
variables,  where  n  is  the  size  of  the  array  and  m  is  a  large 
number  of  closely  spaced  angles  spanning  all  angles.  The  rows 
of  this  matrix  are  passed  through  a  low  pass  filter  to  achieve 
the  0.98  correlation  at  one  beamwidth.  Eor  the  selected 
parameters,  this  limits  the  maximum  achievable  sidelobes. 

In  Eig.  4  we  examine  element  15  failure  in  the  64  element 
uniform  linear  array.  Plotted  are  the  algorithm  results,  the 
original  with  no  failed  elements,  the  original  with  element  15 
set  to  0,  plus  an  optimum  explained  below.  We  used 
beamformer  antenna  patterns  with  the  lowest  possible 
sidelobes  due  to  the  element  pattern  differences  from  element 
to  element  and  thus  useG2(0)  with  procedure  1,  ASlgo^ti  = 
0  dB  (i.e.,  unchanged  sidelobes)  and  =  —3.2  dB.  The 

relevant  parameters  K,  y,  and  taper  loss  are  indicated.  The 
mainlobe  is  a  little  wider  but  the  sidelobes  are  about  the  same 
as  the  pattern  with  no  failed  elements. 

In  order  to  evaluate  the  performance  of  the  correction 
algorithm,  we  include  in  Pig.  4  a  comparison  with  an  optimum 
approach.  We  define  the  optimum  beamformer  as  the  one 
resulting  from  (1)  with  no  manifold  errors  and  y  selected  to 
maintain  unchanged  sidelobes  and  2d  chosen  to  be  the  angular 
width  of  the  K  beams  in  the  algorithm.  Notice  that  the 
algorithm  and  optimum  patterns  virtually  overlay  with 
identical  taper  losses.  The  (All)  bounds  true  sidelobes  at  5.3 
dB  above  the  estimation,  although  the  actual  discrepancy 
between  truth  and  estimation  is  less  than  1  dB. 

Pig.  5  repeats  the  previous  example  except  the  original 
sidelobe  levels  have  been  raised  by  5  dB  by  adding  a  random 
vector  and  thus  we  use  Here  again,  the  algorithm  and 

the  optimum  are  virtually  identical  with  similar  taper  losses. 
The  (A2)  bounds  true  sidelobes  at  4.7  dB  above  the  estimate 
but  again  the  actual  discrepancy  is  less. 

Pig.  6  is  the  same  as  Pig.  4  except  with  three  failed  elements 
15,  32,  53.  Here  again  the  optimum  and  the  algorithm  agree 


well.  The  high  first  sidelobes  are  primarily  due  to  deleting  an 
element  in  the  middle  of  the  array.  Had  we  plotted  results  for 
the  single  element  32  failure,  a  similar  high  first  sidelobe 
appears.  The  bound  given  by  (All)  is  8.3  dB  which  is  not 
very  tight  due  to  the  large  value  of  K. 

Pig.  7  shows  the  amplitudes  of  each  component  of  w(0) 
from  Pigs.  4  and  6  for  both  the  algorithm  and  optimum  which 
are  very  similar.  The  jagged  nature  is  due  to  the  modeling  of 
the  antenna  element  errors.  In  all  cases  maximum  use  is  made 
of  the  largest  contiguous  interval  of  operating  elements. 

Pig.  8  repeats  Pig.  5  except  we  choose  parameters  to  reduce 
the  taper  loss  by  allowing  for  higher  sidelobes.  We  apply 
procedure  2  and  set  TLg^^i  =  2.5  dB  and  =  1  dB. 

The  actual  achieved  change  in  sidelobes  was  1.6  dB.  The 
reduction  in  the  taper  loss  resulted  in  high  first  sidelobes.  Here 
(A2)  bounds  true  sidelobes  at  8.9  dB  above  the  estimate. 

These  examples  demonstrate  the  ability  to  trade  sidelobes 
for  lower  taper  loss. 

B.  Two-dimensional  array  examples 

Here  we  show  the  application  of  the  algorithm  to  a  two- 
dimensional  array.  The  array  errors  were  modeled  in  the  same 
way  as  the  one-dimensional  example  with  errors  correlated  in 
both  dimensions.  Pig.  9  shows  a  low  sidelobe  pattern  for  a  16 
X  16  element  array  having  sidelobe  levels  of  -52  dB  off  the 
cardinal  axis  and  -39  dB  on  the  cardinal  axis.  Pigs.  10  and  11 
show  the  patterns  with  two  failed  elements  at  coordinates  [4,8] 
and  [8,12]  before  and  after  correction  where  the  goal  was  to 
match  the  original  sidelobe  levels.  Here  we  chose  procedure  1, 
and  a  two-dimensional  grid  of  beams  spaced  by  half 
beamwidth.  The  criteria  were  met  with  K  =  13  and  y  =  0.43. 
Pigs.  12  and  13  show  the  beamformer  amplitudes  applied  to 
each  element  before  and  after  correction  where  the  jagged 
nature  is  due  to  the  nature  of  the  antenna  manifold  errors. 
Notice  in  Pig.  13  the  algorithm  places  emphasis  on  the 
elements  in  the  upper  right-hand  corner  to  take  advantage  of 
the  largest  contiguous  area  of  operating  elements.  The 
optimum  approach  does  similarly.  The  achieved  sidelobes  in 
Pig.  1 1  are  unchanged  off  the  cardinal  axis  but  are  raised  2  dB 
on  the  cardinal  axis.  The  taper  loss  has  increased  from  -1.8  dB 
to  -3.0  dB.  Better  taper  losses  are  achievable  with  larger 
values  of  K.  Por  example  with  K  =  25  it  is  possible  to  achieve 
a  taper  loss  of  -1.9  dB,  comparable  sidelobe  levels,  but  higher 
near  in  sidelobes. 

IV.  Sidelobe  Adjustments  with  no  failed  elements 

Por  some  applications  it  may  be  desirable  to  dynamically 
control  the  taper  loss  or  sidelobe  level  with  no  failed  elements. 
Pig.  14  shows  results  using  procedure  2  where  the  taper  loss 
has  been  reduced  by  0.9  dB.  The  pattern  has  high  first 
sidelobes,  near  in  sidelobes  have  been  raised  and  far  out 
sidelobes  are  about  the  same.  Different  values  of  K  and  y  yield 
different  changes  to  the  pattern. 

V.  Pulse-Doppler  radar  application 

As  is  frequently  the  case,  techniques  developed  in  either  the 
spatial  or  the  temporal  domain  frequently  finds  application  in 
the  other  domain.  Consider  a  pulse-Doppler  radar  in  which 
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one  or  more  pulses  are  severely  interfered  with,  but  low 
Doppler  sidelobes  are  needed  with  these  pulses  dropped. 
Mathematically,  the  missing  pulses  replace  the  failed 
elements,  and  the  Doppler  filters  replace  the  low  sidelobe 
beamformers.  The  techniques  in  this  paper  offer  a  rapid  and 
predicable  outcome  for  taper  losses  and  Doppler  sidelobe 
level.  Since  here  manifold  errors  are  likely  not  to  be  an  issue, 
it  may  be  preferable  to  use  (3)  instead  of  the  approximation  in 
(7)  for  the  covariance  matrix.  Eurthermore,  the  techniques  in 
[1-8]  are  also  applicable  but  may  be  more  computationally 
intense  with  little  added  benefit. 

Similarly  temporal  samples  in  the  range  domain  may  be 
interfered  with  and  low  range  sidelobes  are  needed  even  with 
these  dropped  samples.  Mathematically  the  pulse  compression 
filter  replaces  the  low  sidelobe  beamformers. 

VI.  Summary 

We  started  with  the  observation  that  a  low  sidelobe 
beamformer  with  failed  elements  could  be  formed  by  a  linear 
combination  of  a  few  beamformers  around  the  direction  of 
interest  constrained  to  have  a  0  at  the  location  of  the  failed 
elements.  To  choose  the  number  of  beamformers,  K,  and  the 
relative  value  of  modeled  thermal  noise,  y,  we  derived  an 
approximation  for  the  change  in  sidelobe  level  and  the  taper 
loss.  We  have  thus  demonstrated  a  robust  method  for 
calculating  a  low  sidelobe  beamformer  with  failed  elements 
where  we  can  control  the  achieved  sidelobe  levels,  the  taper 
loss,  or  the  width  of  the  mainbeam  region. 

The  approach  is  also  applicable  to  forming  beams  with  no 
failed  elements  to  reduce  taper  losses  by  raising  the  sidelobes. 
Also  discussed  was  the  application  of  the  algorithm  in  pulse- 
Doppler  radars  for  missing  pulses  or  missing  range  samples 
with  the  goal  of  maintaining  low  Doppler  or  range  sidelobes. 

The  method  has  been  shown  to  be  near  optimal  and  thus 
may  find  application  even  for  systems  with  good  knowledge 
of  the  antenna  manifold. 

Appendix 

Here  we  examine  the  two  variations  of  the  sidelobe 
estimator  and  derive  bounds  on  how  much  they  can 
underestimate  the  true  sidelobe  levels. 

Ideally  we  would  like  perfect  agreement  of  estimation  with 
truth.  This  is  impossible  without  detailed  measurements,  but 
for  a  good  choice  of  G^st  we  can  at  least  avoid  the 
unacceptable  situation  of  choosing  a  beamforming  weight 
vector  for  which  SL^stC^)  <<  ^^trueC^)- 

While  the  results  can  be  generalized,  for  simplicity  we 
assume  here  that  sidelobes  of  all  the  beams  Wj^  are  all 
approximately  equal  to  a.  All  w  vectors  in  this  appendix  are 
constructed  from  the  columns  of  the  matrix,  and  so  each 
such  w  uniquely  defines  a  corresponding  c  E  by  the 
equation  w  =  W^c,  or  equivalently,  c  = 

We  define  the  vectors  ej^  by  Wj^  =  Thus  ej^  is  simply 

the  standard  basis  vector  of  that  is,  a  vector  with  a  1  in 
the  position  and  0  elsewhere.  Eor  convenience,  we  let  the 
index  k  vary  from  -  to  +  Since  the  look  direction  of 
interest,  0,  is  fixed,  we  use  the  simplified  notation,  Wq  = 


w(0),  and  Cq  =  (W^Wj^)~^W^Wq.  We  use  the  notation 
Ra  —  ^K'^a  to  represent  the  constraint  vector  in  the  M4: 
beamspace.  When  we  refer  to  taper  loss  in  this  appendix,  we 
always  mean  with  respect  to  the  assumed  steering  vector, 
which  for  the  purposes  of  this  paper  should  be  nearly 
indistinguishable  from  the  true  (unmeasurable)  taper  loss.  We 
remind  the  reader  that  M,  hence  also  G,  is  formed  by 
excluding  a  mainbeam  guard  region  wide  enough  to  exclude 
the  mainbeam  regions  of  all  the  beams  in  W.  Let  G^^t  be  the 
model  covariance  matrix  used  in  sidelobe  estimation,  in  the 
beamspace  basis. 

The  bounds  will  make  use  of  the  following  well-known 
inequality:  If  a^,  a2; are  real  numbers,  then 

(ci^  +  •••  +  <  K  •  +  •••  +  (Al) 

This  inequality  may  be  proven  several  ways.  Eor  example,  one 
can  multiply  out  the  terms  on  the  left-hand  side  of  (Al)  and 
repeatedly  apply  the  simpler  inequality  ZaiUj  <  af  +  of, 

which  follows  easily  from  (a^  —  >  0. 

A.  Diagonal  Model  Covariance  Matrix 

Let  aj^  denote  the  sidelobe  power  levels  of  the  beamformer 
W/^,  and  define  G^^t  =  If  we  are  estimating  the 

sidelobes  of  w  =  and  c  =  then  the  highest 

possible  sidelobes  occur  when  the  terms  in  the  sum  add  in 
phase: 

(^|cj  •  (A2) 

One  can  directly  apply  (A2)  to  compute  a  worst-case  bound  on 
sidelobes  of  a  c  obtained  from  one  of  the  algorithms  in  this 
paper.  In  general,  applying  (Al)  to  the  sum  yields 

Sltruejc)  ^  (Zkfcl  •  ^  (A3) 

^/c 

SO  that,  for  small  values  of  K,  estimation  error  is  acceptably 
low.  Of  course,  estimates  of  average  sidelobe  levels  in  any 
spatial  region  of  interest,  as  well  as  peak  sidelobe  estimates, 
can  be  bounded  in  a  similar  manner. 

In  the  body  of  the  paper  (6)  assumes  the  special  case  in 
which  aj^  =  a  is  independent  of  k.  Examples  using  (6)  apply 
this  bound. 

B.  Model  Covariance  Matrix  For  Wq  Optimal 

Here  we  follow  the  definition  from  Section  II. D.  We  will 
considerably  simplify  our  computations  with  the  following 
additional  notation:  TLj^  =  \qa^k\^  =  l^a^/cl^  is  the  taper 
loss  for  the  k^^  beam,  =  e^Gej^/TLj^,  and  Aj^  i  =  p^  —  pi. 
We  will  sometimes  also  use  \\c\\q  =  c^Gc. 

Using  our  compressed  notation: 

^est  ~  Po  ’  RaRa  "f  tlisg(TL/^  •  (A4) 

It  can  easily  be  verified  that  G^st  has  the  following  desirable 
properties:  1)  G^st  is  a  positive-definite,  Hermitian  matrix,  2) 
ekCest^k  =  efcGgfe  for  all  k,  and  3)  GestSo  =  ia/eo<la)<la- 
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The  first  is  a  requirement  for  a  sensible  model  covariance 
matrix.  The  second  says  that  G^st  agrees  with  G  on  the  basis 
vectors  (known  from  observation).  The  third  says  that 
^0  ^  that  is,  our  original  weight  vector  Cq  can  be 

obtained  by  covariance  matrix  inversion  using  G^^t  in  the 
beamspace,  thus  G^^t  is  a  model  covariance  matrix  for  which 
Wq  is  optimal. 

Now  we  bound  the  error  incurred  by  this  approximation. 
Denote  by  =  Wj^c^pf  the  true  unit-normed  optimal  under 
M  when  restricted  to  the  M4:  beamspace,  so  that  Gc^p^  oc  q^. 

Let  TLopf  =  \qaCopt\^  he  the  taper  loss  of  the  true  optimal, 
and  similarly  define  and  ^k^opt- 

We  will  use  the  following  shortly:  If  x  is  any  vector  such 
that  q^x  =  have 

c^ptGx  =  {GCoptT X  oc  q^x  =  0  (A5) 

Define  vectors  and  8^  by 


^“7^'  (A6) 

Ha^O 

By  these  definitions  one  can  easily  verify  that,  for  all  /c, 

QaTfc  =  0  = 

„  _  (A7) 

Ha^O 

We  will  need  the  following,  which  is  computed  using  (A5), 
(A6),  and  (A7): 

W'^kWc  —  ‘  ^k,opt  (A8) 

Using  (A7)  and  (A8)  we  get  the  following  bound,  which  we 
will  use  shortly: 


('nnt  T  u  T 


—  TLq  •  Popt  T 


T  ^  '  ^k8k 

k^O  G 

||2 

■  ^  ^/c 


<7’L,-p„pt+  ^  l|TollG+^kfe|-||«fc 

\\Ha  0 

<7’Lc-po  +  (^JCfe|-|l«fcllGj  + 
+2(TL,-Ao,opt/^'^|cJ-|l«kllG 


It  remains  to  estimate  the  sidelobes  and  taper  loss  of  w^pf. 
In  the  primary  case  of  interest,  average  sidelobes  cannot  be 
significantly  reduced  from  that  of  Wq  due  to  antenna  errors,  so 
we  take  c^p^GCopt  ~  ocq-  Lor  an  approximation  of  taper  loss, 
note  that  largest  (best)  possible  value  for  T L^p^  is  achieved  by 
orthogonal  projection  of  the  assumed  steering  vector  into 
the  Wj^  subspace.  This  estimation  is  overkill,  but  typically 
adequate.  Substituting  these  estimates  into  our  previous 
definitions  gives  us: 


TL^  v’iw^iw^w^ym^va 

Substituting  (A  12)  and  (A9)  into  (All)  allows  us  to  compute  a 
numerical  bound  on  true  sidelobes,  and  thus,  bound  the 
estimation  error.  In  the  examples  when  (7)  is  used  this  bound 
is  computed. 

We  can  derive  a  much  simpler  bound  in  the  special  limiting 
case  that  «  Po-  Applying  (Al)  to  (All)  gives: 


li  <  Ty■{p,  +  ^,,,p,)  +  KY\cu\^■\\8u\\l 

k^O 

^  rL,-po  +  /^^|cJ2.||5J|2 


—  '^^k  '  (^/c,opt  "f  ^0,opt  "f  ’^\,opt^k,opt) 


Now  we  are  ready  to  bound  -  by  expressing  c  in  the 

c  Gest^ 


{Cq,  Sj^}j^^q  basis.  We  write: 


0  +  ^  Ck^k 


where  =  ^k^  is  simply  the  component  of  c  in  the 
standard  basis.  To  bound  \\c\\q  =  c^Gc  we  will  use  (A5)-(A8) 
and  (A  10),  obtaining 


When  Ag^opt  «  Pg  in  this  case  it  is  typical  to  also  have 
^o,opt  <<  ^k,opt  for  k  ^  0,  allowing  us  to  approximate  (A9)  as 
||5fc|li  ^  rLfc  -Afc.o,  yielding 


c«Gc  <  TL,  •  Po  +  •  TL„  ■  A, 


In  comparison,  our  estimate  (A4)  gives 
C^GestC  =  TLc  -Po  +^kfck  -T 


y  ■  ^k.o  (A15) 


Therefore  in  this  case  the  bound  may  be  approximated  simply 
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c^Gc 

c^Gp^tC 


<  1  +  (K  -  1) 


where  (Tq  =  —  and 

«o 


TL{c) 


(A16) 
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n  ^  2(0) 

antenna  ^ 

elements  T 

w(0) '  '  beamformer 

Fig.  1.  Beamformer  block  diagram. 


Fig.  2.  Procedure  for  setting  a  goal  for  sidelobes  and  bound 
on  taper  loss. 


Fig.  3.  Procedure  for  setting  a  goal  for  taper  loss  and  a 
bound  on  sidelobes. 
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Fig.  4.  Antenna  patterns  with  failed  element  15  with  62(0). 


sin(angle) 

Fig.  5.  Antenna  patterns  with  failed  element  15  with  Gi(0). 


sin(angle) 


Fig.  6.  Antenna  patterns  with  failed  elements  15,  32,  and  53 
with  62(0). 


10  20  30  40  50  60 

Element  # 

Fig.  7.  Weight  vector  amplitudes  for  cases  in  Figs.  5  and  7. 
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sin(angle) 

Eig.  8.  Antenna  patterns  with  taper  loss  limited  to  -2.5  dB. 
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Eig.  9.  Low  sidelobe  pattern  for  16  x  16  array. 
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Eig.  1 1 .  Corrected  pattern  with  two  failed  elements. 
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Eig.  13.  Beamformer  amplitudes  with  failed  elements  [8,4] 
and  [12,8]. 
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sin(angle) 


Eig.  14.  Antenna  patterns  with  no  failed  elements  and 
reduced  taper  loss. 


